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1. Introduction 

In its simplest manifestation, gravitomagnetism is caused by the current of mass- 
energy in direct analogy with electrodynamics. In particular, gravitomagnetic effects 
can be found in the external gravitational field of rotating astronomical masses. 

In this article we are interested in the exterior Kerr-Newman-Taub-NUT (KNTN) 
solution given in the Boyer-Lindquist coordinates by the metric 1,2,3,4,5 

1 2 

ds 2 = - - (A - a 2 sin 2 9)dt 2 + - [Ax - a(S + ax) sin 2 6»]did0 

+ -[(£ + a X ) 2 sin 2 9 - x 2 A]d0 2 + ^dr 2 + £d<9 2 (1.1) 
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and the corresponding electromagnetic Faraday tensor can be expressed in terms of 
the 2-form 

F = $^{[r 2 - {i + a cos 9) 2 }dr A (dt - X d(j)) 

+ 2r sin 9{i + a cos 9)d6 A [(r 2 + a 2 + £ 2 )d<j) - adt]}. (1.2) 
Here £, A, and x are defined by 

£ = r 2 + (£ + acos9) 2 , A = r 2 - 2Mr - I 2 + a 2 + Q 2 , x = asm 2 9 - 2£cos9 . (1.3) 

Units are chosen such that G = c = 1, so that (M,Q,a,£) all have the dimension 
of length: the source has mass M, electric charge Q, angular momentum J = Ma 
(i.e., gravitomagnetic dipole moment) along the z-direction, and gravitomagnetic 
monopole moment n = —£, where £ is the NUT parameter. In this family the charged 
KNTN metric can be obtained from the uncharged Kerr-Taub-NUT (KTN) metric 
by the simple transformation M — > M — Q 2 /(2r) in A, while the Taub-NUT metric 
corresponds to a = 6 7 8 9. It follows from equation H1.2|l that in the far zone 
(r 3> a and r ^ \£\), the lowest-order contribution to the electric field is the radial 
Coulomb field, while the lowest-order contribution to the magnetic field is due to the 
field generated by the magnetic dipole moment Qa of the source as well as terms of 
similar order proportional to 2Q£. 

An interesting example of a gravitomagnetic phenomenon is provided by the clock 
effect for circular equatorial orbits around rotating masses |E3[n]. To illustrate the 
situation in a more general context, consider the circular equatorial motion of a test 
particle of mass m and negative electric charge — q < around the Kerr-Newman black 
hole of mass M, charge Q > 0, and angular momentum J = Ma. Here we neglect 
all (electromagnetic and gravitational) radiative effects for the sake of simplicity. A 
Newtonian analysis implies that the test particle moves along a circle of radius r with 
uniform angular frequency given by (^-pr^) 1 ^ 2 , where —rj = —q/m < is the charge- 
to-mass ratio of the test mass. This Newtonian result is expected to emerge from the 
following general relativistic analysis when the test particle is very far from the black 
hole. In the relativistic theory, the Lorentz equation of motion 

dx a dxP „,. dx u 



must be solved. Here r is the proper time of the test particle along its world line. A 
detailed analysis reveals that in this case Eq. I|1.4I) reduces to 

where ±u>k (u>k > 0) defined by 



M 



y2 



is the Keplerian frequency for a neutral test particle in circular orbit around a Reissner- 
Nordstrom black hole and ±wc (wc > 0) defined by 

ul = r)% (1.8) 
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would be the frequency of motion of the test particle of charge —r\m if only the 
Coulomb interaction were present. 

To analyze the system (I1.5fl - <|1.6[1 . first set a = 0; then (dtp / dt) 2 satisfies a 
quadratic equation that has two possible solutions. These can be distinguished by 
the fact that in one case dt/dr > and in the other case dt/dr < 0. We must impose 
the physical requirement that the time coordinate t increases along timelike world 
lines. With this postulate, one then obtains d<fi/dt = ±fl, where fl > is given by 

n 2 = uj 2 k - ^r 2 oj^, + Nlu 2 c (1.9) 



with 



V r r z 8 

Here VI 2 — in contrast to the second solution that we have ruled out — properly reduces 
to the Newtonian result far from the black hole. 

Next consider a charged rotating black hole (a ^ 0,Q ^ 0). The results can be 
simply stated for a/M <C 1, in which case to first order in a/M one finds 

dt , 1 , 

d<r^ +a(1 - s) ' (L11) 

where 

Q=l VQ. (112) 

2 N [M-91(i + l 1 f ) + v Q N ] ■ 

Thus if t + denotes the coordinate time period for a complete revolution around 
the black hole in the counterclockwise (clockwise) sense (A</> = ±27r), then 

t + -t-= 4%a(l - Q) , (1.13) 

which is the expression for the single-clock clock effect [11)11111 in this case. Note that 
as r — > oo this result reduces to 

M + ±r)Q , 

*+-<- = 4 ™mTW' ^ 

so that the effect is always less than Aira, the result for a neutral particle. It is 
interesting to investigate the two-clock clock effect ^3^3 m this case as well, i.e. 
t + — r_, where r is the corresponding proper time as measured along the geodesies. 
One finds 



dr . n 



2 



= ±^L + < 1-Q>), (1.15) 



where 



It follows that 



Q'= 1 2 N ^( 1 + ^)- (1-16) 



r+ - t_ = 47ra(l - Q!) (1.17) 

and that Q' — > Q as r — > oo. These results, valid only to first order in a/M, indicate a 
certain gravitomagnetic temporal structure around the Kerr-Newman black hole that 
is consistent with the axial symmetry of the underlying spacetime structure. It is the 
purpose of this paper to clarify how gravitomagnetism is affected by the presence of 
the gravitomagnetic monopole moment [i = —£, which has spherical symmetry. 
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In the following section, the gravitomagnetic field is studied and the results are 
applied in section 3 to the propagation of test electromagnetic waves in the Kerr-Taub- 
NUT spacetime. Section 4 describes circular holonomy in this spacetime. Spherical 
orbits are considered in section 5. This is followed in the final section by a discussion 
of our results. 



2. Gravitomagnetic field 

Even though the spacetime curvature vanishes asymptotically at large spatial 
distances, the KNTN spacetime is not asymptotically flat when £ ^ 0, a fact intimately 
related to the gravitomagnetic monopole character of the source associated with this 
parameter. At present there is no evidence for the existence of either magnetic or 
gravitomagnetic monopoles. Moreover, the lack of asymptotic flatness for a spacetime 
describing a localized source appears to be unphysical. However, like the Godel 
spacetime, the Kerr-Newman-Taub-NUT spacetime can play a useful role in clarifying 
important gravitomagnetic features of classical general relativity. 

With this goal in mind, we set Q = and expand the metric to first order in a 
and p = in this way one can more easily interpret both parameters in the context 
of gravitomagnetism. One finds 

ds 2 = - (l - —) dt 2 + ( 1 - —) 1 dr 2 + r 2 (d9 2 + sin 2 9 dcj?) 



aM ■ 2 a a A 2M \ 
sin 6 - fj, cos 8 I 1 I 



dtd<j> , (2.1) 



namely, the Schwarzschild metric plus the Lense-Thirring and Taub-NUT terms. 
While the time-coordinate slices are intrinsically asymptotically flat, the fact that 
g tl f, — ► — 2£cos0 as r — > oo implies that the spacetime is not asymptotically flat [HJ[7]- 
By introducing the isotropic radial coordinate p defined by r = p(l + ( f>/2) 2 , 
where $ = M/p is the (sign-reversed) Newtonian gravitational potential, and the 
corresponding isotropic Cartesian-like coordinates x — p sin 9 cos 4>, y — p sin 9 sin <j), 
and z = pcos9 (in 3-vector form p = pp = [x, y, z)), the metric takes the form 



ds = ~ {iT$/2 ) dt + $/2) {dx + dy + dz ) 

- 4A (g ) • dpdt . (2.2) 

We interpret this metric as representing a background Minkowski spacetime in 
inertial coordinates (t, x, y, z) together with a source characterized by a gravitoelectric 
potential $ and a linearized gravitomagnetic vector potential A( g ) given by 

A( g) = T(J x p) <-> A (g ) • dp = Tp 2 sin 2 9 d<p (2.3) 

(i.e., the slicing gravitomagnetic potential 1-form N a dx a = g ta dx a when multiplied by 
—2, see 52])- Here J = Mai is the angular momentum, 



aM pz /l-$/2 



2 



T p 3 (l + $/2) 2 p(x 2 + y 2 ) \1 + <P/2J ' ^'^ 
and the 3-vector operations are carried out in the flat spatial metric, so that 

(z X p) ■ dp = p 2 sin 2 9 d<j) = (p 2 - z 2 )d<\> . (2.5) 
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In the case that M/p <C 1 far from the source, one has 

Ma fxz J x p pzzx p 

~ p 3 p(.T 2 +j/ 2 )' (g) ~ p 3 p(2; 2 +?; 2 )' 1 j 

and the gravitomagnetic vector field B( g ) = V x A( g ) is given by 

B (g) = ^[3(p.J)p-J] + ^. (2.7) 

Thus J corresponds to the gravitomagnetic dipole moment and p to the 
gravitomagnetic monopolc moment. The threading gravitomagnetic vector potential 
(i.e., the 1-form M a dx a = — (gt a /gtt)dx a , see J2]) defines another vector 

^ A + $/2\ 2 - „ Jxp „ iizzx p . 

7 = " 2 ^(T^j Jx ^- %(p-m/2) 2 +2 ^T7)' ^ 

so that 7 — 2A( g ) for p ^> M. 

One can now use the metric (I2.2|l to treat phenomena in terms of a gravitoelectric 
field E( g ) = — V$ and a gravitomagnetic field B/ g ) = — V*, where 

* = ^ + ^ (2.9) 

is the gravitomagnetic scalar potential I n the definition of gravitoelectric 

and gravitomagnetic fields, certain conventional numerical factors of 2, etc., are 
unavoidable depending on the specifics of the analogy with electromagnetism. It 
turns out that in the description of the linearized gravitational motion of a test 
particle of inertial mass m around a central source of inertial mass M, one can 
employ the same formulas as in classical electrodynamics if one assigns a gravitoelectric 
charge M and a gravitomagnetic charge 2M to the central source and corresponding 
gravitational charges of —m and —2m to the test mass ^31- The different signs 
ensure that gravitation is attractive; moreover, the ratio of gravitomagnetic charge 
to the gravitoelectric charge is always 2 due to the spin-2 character of the linearized 
gravitational field. We employ this convention in the present work. 

Geodesic motion of test particles in the field of a gravitomagnetic monopole has 
been considered in |14U15| and in references cited there, while lensing aspects of a 
gravitomagnetic monopole have been discussed briefly in |15| . In the following section 
we consider the effect of this monopole on the propagation of electromagnetic waves. 



3. Wave propagation 

Consider the propagation of electromagnetic waves in the Kerr-Taub-NUT spacetime. 
We work to linear order in a and p and express the electromagnetic perturbations 
in terms of the background Cartesian-like coordinates (t, x, y, z) introduced above. 
Writing the sourcefree Maxwell equations in these coordinates as f[ Q /3, 7 ] = and 
[y/— gF al3 ].f3 — and introducing the standard decompositions F a p — > (E,B) and 
^J—gF a P — > (— D,H) (TO], one finds the equivalent flat spacetime Maxwell equations 
in inertial coordinates in the presence of an effective material medium [161117) . i.e., 

V • D = , V • B = , VxE = -d t B , V x H = 9 t D . (3.1) 

Therefore we may think of the electromagnetic field as propagating in an inertial frame 
in Minkowski spacetime but in the presence of a gravitational "medium" with certain 
constituitive properties given by |16U17lll8lll9j 

A = tijEj - [7 x H]i , B, = p l3 H 3 + [7 x E]< , (3.2) 
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Figure 1. Propagation of null rays along a closed contour around a mass M 
with gravitomagnetic monopole (fi = —£) and dipolc (J = Ma) moments. The 
"mirrors" m\-m/± in this schematic diagram could be transponders on-board 
artificial satellites. 



where e,j = /iy = — \/ ~ gg 11 1 gtt and 7^ = —gu/gtt ( = Mi, the threading shift 1-form 
as above). For the present metric one has 

ur K r (l + ^/2) 3 .... 

<Hj = fiij = AfSij , N = 1 _ ^ 2 , (3.3) 

and 7 is given by Eq. Q2.8fl . Here A/" can be interpreted as an effective index of 
refraction of the gravitational medium and 7 represents the gyrotropic aspects of the 
medium. We note that 7 is singular along the entire z-axis due to the presence of the 
term involving fi 0. However, this singularity can be removed by a local coordinate 
transformation in any given exterior neighborhood , a fact that is consistent with 
the gauge dependence of the vector potential. 

Since these constitutive quantities are all independent of time, the electromagnetic 
waves propagate with constant frequency u> in this gravitational medium. With a 
time dependence of the form exp(— iuit), Maxwell's equations then reduce to the wave 
equation 

(i _1 V -w7)xW = -iwAfW (3.4) 

for the complex Kramers vector W = E + iH. Note that the gravitomagnetic 
contribution to this Dirac-type equation enters via the substitution i _1 V — > i _1 V— ujj. 
Therefore, it follows from the Aharonov-Bohm effect that radiation following a closed 
spatial path C — via reflection from devices such as mirrors, for example (see Fig. 1) — 
in the exterior field of a gravitomagnetic source exhibits the following phase difference 
comparing the two waves after one loop in the counterclockwise (+) direction and one 
loop in the clockwise (— ) direction along the closed path C ^21 



y + -V.=2uf 7 • dp . (3.5) 

For p M, V x 7 ~ — 2B(g); therefore, one finds that <p+ — cp- is —Alo times the 
gravitomagnetic flux threaded by the closed path C. 

If we consider matter waves propagating around this path, corresponding to 
radiation of particles with nonzero rest mass m > instead of photons, then to in 
Eq. I|3.5|l must be interpreted as the de Broglie frequency of the matter waves, i.e. 
lo — > m/h. Moreover, a different treatment is required if instead of C a path defined 
by optical fibers (or a ring laser) is employed. Clearly the constitutive properties of 
the corresponding medium should then be taken into account as well. 

If the dominant wavelength of the waves is negligible compared to the 
characteristic length scales of the medium, one can interpret Eq. I|3.5|l in terms of 
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rays of radiation. In this case, let P be an observer at a point on the closed path C. 
Then according to this observer the interval of the coordinate time t + (t_) for the rays 
to go around C in the positive (negative) sense is such that ip + — (p_ = uj{t + — 
Thus 



t+-t-=2j> 7- dp , (3.6) 

which is reminiscent of the Sagnac effect via the gravitational Larmor theorem |13| . 

Consider now a closed path C in the (x, y)-plane far from the source (p ^> M) 
and note that the monopole contribution to t + — t_ in Eq. (|3.6|) vanishes due to the 
fact that z = in this plane (see Eq. (|2.8|0 . It should also vanish in any other plane 
through the origin by the spherical symmetry of the monopole. This is clear from a 
simple application of Stokes' theorem (see below) if C does not encircle the singularity 
at p = 0. If it does, then let Co be a circle in this plane around the origin, such that 
C is always outside Co. It follows from Stokes' theorem that 

z[zxp].d P = rp_ds =o 

c-c x z + y J P A 

since p in this plane is orthogonal to the area element dS by definition. It now remains 
to calculate 

J Co x + y 

in this plane and to show that it vanishes regardless of the radius r of the circle Co . 
Taking advantage of the axial symmetry of the "medium" , it suffices to calculate / in 
the (x', y')-plane, where the (x 1 , y' , z 1 ) coordinate system is obtained from the (x, y, z) 
system by a rotation of angle a about the x-axis. Thus 

x = x, y ~ y' cos a — z' sin a, z = y' sin a + z' cos a. (3-9) 

A straightforward calculation then results in 

C sin0d</> „ ,„ 

1 = smacosa / ^ ^ — = " , (3.10) 

Jo cos 2 a + sin a cos 2 <f> 

where x' — ro cos^>, y' — sin^>, and z' = along Co- 

It is interesting to point out that an equation very similar to Eq. H3.4JI would hold 
in a background with Cartesian-like coordinates for the propagation of electromagnetic 
waves in the general Kerr-Taub-NUT spacetime. In this spacetime, one must identify 
t with t + n(8ir£ ) for any integer n; that is, t is assumed to be periodic with period 
8ir£ P3EJ- Thus in the presence of a gravitomagnetic monopole moment a temporal 
delay equation of the form (3.6) should be interpreted modulo an integer multiple of 
8ir£; moreover, a time dependence of the form exp(— icot) implies that 4uj£ must be an 
integer. Hence the photon energy is quantized in units of Ti/(A£). Various aspects of 
this quantization condition have been explored in |5ni|2]. We are interested in very 
high frequency radiation; therefore, this discreteness will not be pursued further here. 

Let us note that in the JWKB limit, one can solve the wave equation l|3.4(l with 
the ansatz 

W = W e -M*-*o)+«^(x) 5 (3.n) 
where Wo = W(to,Xo) and 

5(x) = f 7 (x') • dx' + f AA(x') k • dx'. (3.12) 

J xo J xo 
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Here k is a constant unit vector that indicates the direction of wave propagation in 
the background inertial frame. It is useful to write the eikonal S in the form 

5(x) = k • (x - xo) + f [jV(x') - 1] k • dx' + f 7(x')-dx', (3.13) 

where the gravitational delay effects are separated from the simple plane-wave 
propagation aspects of the solution. Moreover, let k, ni, and A2 form an orthonormal 
triad in the background Euclidean space such that k x hi = n 2 . Then W = 
A(hi + m 2 ), where A is a complex amplitude. 

In terms of ray propagation, we may interpret these results as implying that a 
ray starting from Xo at to will reach x at t such that 

t-to = \x-x \+5t GE + 5t GM , (3.14) 

where StcE is the gravitoelectric time delay given by 

rV(x')-l]k-dx'. (3.15) 

For $ << 1, TV ^ 1 + 2$ and St GE reduces to the integral of 2$ along the straight line 
from Xo to x, which is the Shapiro time delay. Moreover, 5t G M is the gravitomagnetic 
time delay given by 

6t GM = f 7(x')'dx', (3.16) 

which vanishes if x — x is in the same plane as the z-axis. We note that the monopole 
contribution to Eq. (|3.16|) . modulo an integer multiple of 8tt£, vanishes in the equatorial 
plane. 

It has recently been shown that the gravitomagnetic time delay caused by 
rotational motion could make a significant contribution to the gravitational lensing 
delay of extragalactic sources and should therefore be taken into account in the 
interpretation of observational data 22, 23, 23- The time delay around a closed loop 
would involve a gravitoelectric as well as a gravitomagnetic component. If the loop is 
traversed in the opposite sense, the gravitoelectric time delay for the stationary field 
under consideration remains the same but the gravitomagnetic time delay changes 
sign. In this way, one recovers Eq. I|3.6|l for the total time delay t + — t_ from a 
different standpoint. 

Finally, let us consider the influence of the gravitomagnetic monopole on the 
polarization of light. The static background Schwarzschild spacetime has no influence 
on the polarization of test electromagnetic radiation due to its spherical symmetry |19| . 
We therefore assume that $ -c 1 and consider linear perturbations of Minkowski 
spacetime such that g^ = rj^ v + h^, 

/too = 2$, /',., 2<!.<v h 0i = 2A\ g) . (3.17) 

Far from the source, V • Am — and □ A( g ) = 0, so that the trace-reversed 
perturbations 

hpu = hfxv - ^hr}^ (3.18) 

satisfy W v „ = and Oh^ = 0. It then follows from the analysis in |19ll25| that in 
a gravitomagnetic field B( g ), the plane of linear polarization rotates at a rate given 

by puna 

§=B (g) -k, (3.19) 
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where C is the rotation angle and k is the unit propagation vector of the 
electromagnetic radiation. This Skrotskii effect JS| is the gravitational analog of the 
Faraday effect in electrodynamics. Thus for electromagnetic radiation propagating 
from p 1 to p 2 , the plane of polarization rotates by the angle ( = ^f(Pi) — ^{P2) m this 
linear approximation. In particular, for a ray propagating from a radius po 3> M out 
to infinity, the plane of the polarization rotates by an angle 

Co = *o = + ± ■ (3.20) 

Po Po 

This result is consistent with the spherical symmetry of the monopole. Moreover, it 
follows from Eq. (|3.20|) that the net angle of rotation of the plane of polarization is 
zero in the scattering case of radiation propagating from — oo to oo. 

To go beyond the JWKB limit described above, one must consider the 
gravitomagnetic phenomenon of helicity-rotation coupling as well; see, for instance, 
|19H2fi| for the case of the Kerr spacetime. 

In connection with the spin-rotation coupling, let us consider the motion of a 
classical spinning test particle in the linearized KTN spacetime given by Eq. (|2.2|) . In 
the absence of the gravitational source, we may assume that the free test particle is 
at rest and carries a constant spin vector with its direction fixed in the inertial frame. 
In the presence of the gravitational source, however, the spin of the particle couples 
to the spacetime curvature resulting in the Mathisson-Papapetrou force 

F a =-^R af 3 J su f3 S^ , (3.21) 

where u a is the 4-velocity and S a ° is the spin tensor of the particle. For a "point" 
particle with spin, we need to impose the Pirani supplementary condition S a ^Uf3 = 
for the sake of completeness. The spin vector of the particle is defined by 

S a = -\{-g) l ^e ot0l5 u S^ . (3.22) 

In the linear approximation with u a (1,0,0,0), S° l « 0, and S 1 ^ « e l; > k Sk, 
we find that in Eq. I|3.21|) the Mathisson-Papapetrou force is Fo — and F,; = 
— (B( g ))j,i5 J . This force may be expressed as the negative gradient of a gravitational 
spin potential energy 7i = S - B( g ) . It turns out that this classical result corresponds in 
the JWKB approximation to wave-mechanical (quantum) results. Using Eq. 1)2.7(1 . we 
see that the energy Ti. consists of two terms: the familiar gravitational spin-rotation 
coupling term (~ 10~ 28 eV for a spin-i particle in an Earth-based laboratory) and 
a new term /xS • pj p 3 that is due to the coupling of spin with the gravitomagnetic 
monopole moment of the source. 

4. Circular holonomy 

Turning now to the parallel transport around closed 0-loops, a key result of |57] for the 
Kerr spacetime involves the relation between the gravitomagnetic temporal structure 
around a rotating mass and circular equatorial holonomy, namely, the holonomy of 
the time component X 1 = dt(X) of an arbitrary vector field X does not in general 
vanish for a spacelike circular equatorial orbit around a rotating mass. This temporal 
holonomy involves a Lorentz boost and the corresponding time dilation indicates the 
existence of the gravitomagnetic clock effect. The clock effect and the temporal 
holonomy vanish in the Schwarzschild spacetime. In the case of the Taub-NUT 
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gravitational field, the corresponding equation of parallel transport around a spacelike 
circular equatorial orbit implies that 

^ + ,*° = 0, (4.1) 

so the temporal holonomy is in general nonzero for £ ^ 0. This indicates that 
the gravitomagnetic monopolc moment of the source produces a separate temporal 
structure around the Taub-NUT source as indicated, for instance, by the Skrotskii 
effect discussed above. 

The gravitomagnetic clock effect is essentially absent in this case due to the 
spherical symmetry of the monopole source, which is a gravitational dyon consisting 
of a gravitoelectric monopole (represented by the mass M) and a gravitomagnetic 
monopole (represented by the NUT parameter £). Moreover, timelike circular 
equatorial geodesic orbits do not exist in the Taub-NUT spacetime. That is, £ 
breaks the degeneracy of the Schwarzschild circular geodesic orbits by lifting them 
off the equatorial plane in opposite directions. This could be intuitively understood 
by considering the monopole contribution to the gravitomagnetic Lorentz force per 
unit mass of the test particle, i.e., 2/ip x v/p 3 28 j. To first order in £, however, 
the orbital frequencies are the same as the Keplerian frequencies in the Schwarzschild 
spacetime, dcj)/dt = ±W( K ) with cj( K ) = (Af/r 3 ) 1 / 2 . This is treated in more detail in 
the following section. 



5. Geodesies in the KNTN spacetime 



The motion of a free uncharged test particle in the Kerr-Newman- Taub-NUT 
spacetime is given by 



1 

± - 
S 



R 



R = -A(m 



2 2 

r 



= ± 
1 

1 

S 



sin 



t = - 



(L z - E X ) - 
a 2 + £ 2 ) K- 



<d = K-m\£- 
P 1 \ 



- K) + P 2 , P = E(r 2 +a 2 + £ 2 ) - aL z 
acos9) 2 - — L-(L,- E X f , 



sin 



X 



sin 



(L z - E X ) 



(5.1) 
(5.2) 
(5.3) 
(5.4) 



where m, E, L z , and K, are the four constants of motion for the particle and represent 
respectively its rest mass, energy, angular momentum about the z-axis, and Carter's 
constant |2H]. Here the dot represents differentiation with respect to the affine 
parameter A whose differential is related to the proper time by dr = m dX when 
r7i 7^ 0. For the timelike orbit of a test particle with nonzero rest mass m, the principle 
of equivalence implies that the orbit is uniquely characterized by the reduced constants 
E = E/m, L z = L z /m, and K. = JC/m 2 . A detailed examination of the solutions of 
these equations is beyond the scope of this article; therefore, we shall only consider 
certain special orbits here. 

Let us first assume that a = 0. In this case, we choose L z ^ and K. = 
(m 2 — 4E 2 )£ 2 + L 2 ] one can show from Eq. (5.2) that 6 is then fixed and the 
orbit is confined to a cone with the opening angle 9 given by cos 9 = —2£E/L Z . 
It follows that in this case the equations of motion on the cone depend on £ only via 
£ 2 . As £ — > 0, such an orbit reduces to an equatorial plane orbit; in particular, for 
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constant r we recover circular equatorial orbits. The reflection of the orbit on the cone 
about the equatorial plane can be achieved by simply changing the sign of L z since 
dc/)/dr = L z /(r 2 + £ 2 ); that is, the sense of motion is reversed below the equatorial 
plane. As the gravitomagnetic monopole moment vanishes, the orbits on the cones 
with cos 9 = ±\2£E/ L z \ degenerate to a single orbit on the equatorial plane. 

The absence of timelike circular equatorial geodesies extends to the KTN 
spacetime. More generally, we can look for geodesic orbits that have constant r and 
9 coordinates. These would be special circular geodesic orbits around the rotation 
axis that for 9 — tt/2 coincide with equatorial circular geodesies. A detailed analytic 
investigation of the geodesic equation shows that for o^O and £ ^ no such (timelike, 
null, or spacelike) geodesies exist in the KTN spacetime. 

Let us next consider timelike spherical geodesies in the KTN spacetime. While 
such stable orbits exist all the way down to the horizon (for the corotating orbits with 
a = M) in the Kerr case [SUH^, they are bounded away from the horizon of the 
Kerr-Taub-NUT spacetime, as can be seen by generalizing Fig. 3 of Wilkins [3U| and 
his related calculations. To first order in £, however, the radius of a spherical orbit 
coincides with that of a Wilkins orbit with the same orbital parameters (E, L Z ,}C). In 
this case the orbit follows a spiraling spherical path that stays within latitudes given 
to first order in £ by 9 ni i n = 9 + 9\ and # max = tt — 9 + 6\. Here 9 £ (0, tt/2) is the 
unperturbed latitude satisfying 

(L z - aE sin 2 9 f = [K - a 2 cos 2 O ) sin 2 9 . (5.5) 

If 9q is a solution of this equation, then so is it — 9q; therefore, the Wilkins orbits 
are symmetric in latitude about the equatorial plane. The first order perturbation in 
latitude 9\ is obtained from 9 = and can be expressed as 

e _ ( t \ 2EL z + a(l-2E 2 )sm 2 9 

1 \sm9 ) JC + 2aEL z -a 2 + 2a 2 (l-E 2 )sin 2 9 Q ' 

Thus to linear order in the strength of the gravitomagnetic monopole moment, the 
spherical orbits are shifted up or down and are no longer symmetric in latitude about 
the equatorial plane. 

6. Discussion 

Following up on our previous work ^U|, we have considered some of the main 
effects associated with the existence of a gravitomagnetic monopole moment. The 
time coordinate t is periodic in the KNTN spacetime with a period that is 
simply proportional to the gravitomagnetic monopole moment of the source. The 
gravitational time delay and the rotation of the plane of polarization of electromagnetic 
waves studied in this paper further elucidate the special temporal structure in 
this spacetime. We have also briefly studied the influence of the gravitomagnetic 
monopole moment on the motion of a spinning particle, and on the spherical 
orbits in KTN spacetime generalizing some of the previous results of Wilkins \','>()\ . 
There is no observational evidence at present for the existence of a gravitomagnetic 
monopole |14II15| . Our work could in principle be combined with the analysis 
of astronomical data in order to set upper limits on the possible existence of 
gravitomagnetic monopoles. 
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